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Abstract 


A novel methodology for the on-line diagnosis of PEMFC, intended to be used in portable control systems, is proposed. First, an equivalent-circuit 
model for PEMFC is proposed. This model has a low computational cost and it allows reproducing the dynamic behavior of the PEMFC. The 
relationship among the parameters of this simplified model and PEMFC electrochemical parameters is stated. Secondly, a procedure to estimate 
all the parameters of this model from the cell’s transient response after current interruption (CI) is proposed. The equipment required to make 
these CI measurements is easily portable and inexpensive. Thirdly, the procedure to calculate some relevant PEMFC electrochemical parameters 
from the equivalent-circuit parameters is discussed. These PEMFC parameters are the double layer capacitance, the diffusion resistance, the charge 
transfer resistance, the diffusion-related time constant, and the membrane resistance. Finally, this PEMFC assessment methodology is applied to 
the diagnosis of the cathode flooding. The diffusion resistance is found to have a dominant effect on the amount of water accumulated inside 
the cell. The described methodology can be used to characterize other phenomena, including the cathode and membrane drying, the membrane 


degradation, and the anode poisoning by carbon monoxide. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Our society is suffering from fossil fuel shortage. Fossil fuels 
(i.e., coal, oil, and natural gas) also contribute to a number of 
environmental problems during their extraction, transportation, 
and use. In particular, their use to obtain energy generates emis- 
sions that contribute to global warming and climate change. As 
an alternative, fuel cells are one of the most promising means of 
producing energy in portable systems. 

A considerable research effort has been made in order to 
understand the physical and chemical phenomena involved in 
the PEMFC operation. Numerous mathematical models have 
been proposed, including 1D [1-4], 2D [5], and 3D [6,7] static 
models; and 1D [8,9], and 3D [10] dynamic models. The 
water flooding in the PEMFC cathode has been modelled from 
physical—chemical first principles [11]. These detailed models 
contain a large number of adjustable parameters and their com- 
putation time is large. As a consequence, these models are not 
well suited for control applications. 


* Corresponding author. Tel.: +34 913987147; fax: +34 913987690. 
E-mail addresses: marubio @dia.uned.es (M.A. Rubio), 
aurquia @dia.uned.es (A. Urquia), sdormido @dia.uned.es (S. Dormido). 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
doi: 10.1016/j.jpowsour.2007.06.072 


PEMFC models with low computational cost, intended for 
use in control applications, have also been developed. Randles 
electric models are frequently used in this context. The combined 
use of EIS and Randles models for over-potential calculation 
was proposed by Wagner [12]. The anode poisoning by carbon 
monoxide was discussed by Ciureanu and Wang [13]. The rela- 
tionship between the amount of water and the cell resistance 
was studied by Andreaus et al. [14], and Yuan et al. [15]. The 
effect of the cathode flooding and drying on the cell behavior 
was addressed in [16]. 

In this context, EIS techniques are employed to study the 
fuel cell behavior, and to estimate the value of the cell model 
parameters [17—20]. However, EIS equipment is too expensive 
and bulky for use on in-field assessment of operating commercial 
cells. As a consequence, this approach to PEMFC diagnosis is 
not suited for in-field control applications. 

An alternative methodology for the in-field diagnosis of 
PEMEC, that can be implemented in portable control systems, 
is proposed in this manuscript. To achieve this goal, a Randles 
electric model of the PEMFC is proposed. This dynamic model 
has a low computational cost and it has been specially formu- 
lated to facilitate the estimation of their parameters from the 
cell response after current interruption (CI). The advantage of 
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Nomenclature 

Ca double layer capacitance (F) 

Ce number of electrons in 1 C (CT!) 

Cg concentration in cathode active layer (mol m7?) 
D diffusion coefficient (m? s7!) 

Ew open-circuit voltage (V) 

F Faraday constant (F) 

j imaginary unit 

My,0 water molar weight (18.0157 g mol!) 
n number of electrons 

Na Avogadro constant (6.0221415 mol`!) 
R perfect gas constant (J mol~! K7!) 

Ra diffusion resistance (2) 

Rm membrane resistance ({2) 

Rp charge transfer resistance (Q) 

s laplace transform variable 

S active area (m7) 

Sa electrode surface (cm?) 

T temperature (K) 

Zw Warburg 

Greek letters 

ô diffusion layer width (m) 

Na anode over-voltage (V) 

Ne cathode over-voltage (V) 

Nm membrane over-voltage (V) 

Td diffusion-related time constant (s) 
xXH,0 water load (g cm~?) 

w frequency (rad s7!) 


this approach is that the equipment required to perform the CI 
experiments is inexpensive and easily portable. 

The relationship among the parameters of this Randles model 
and the following PEMFC electrochemical parameters is dis- 
cussed: the double layer capacitance, the diffusion resistance, 
the charge transfer resistance, the diffusion-related time con- 
stant, and the membrane resistance. The procedure to estimate 
the value of these physical parameters from the Randles model 
parameters is proposed. 

The use of CI techniques for cell state characterization is not 
new. Some cell diagnosis methods, based on the CI technique, 
have been proposed by other authors. In particular, CI techniques 
can be successfully applied to the estimation of the FC resistance 
[21]. Other methods, based on the analysis of the cell response 
after a step change in the load, are used to validate the data 
obtained from the EIS study [22]. 

As a difference with other authors’ methods, the proposed 
methodology allows estimating the significant operation param- 
eters of the PEMFC only from the CI data. Rich information 
about the cell dynamic behavior is contained in the cell response 
to step and impulse inputs. In this case, the information is 
extracted from the cell response after CI. The impulse has a wide 
spectrum of frequencies, exciting the FC fundamental phenom- 
ena [23-25]. 


2. Experimental set-up and procedure 


In order to study the effect of the cathode flooding on the 
cell model parameters, the parameter estimation methodology 
proposed in Section 3 is repeatedly applied during the flooding 
process. The data obtained from each CI experiment run, corre- 
sponding to a given cathode flooding level, is used to estimate 
the parameters of the cell electric model. The experiment results 
are discussed in Section 4. 


2.1. Experimental set-up 


PEMFC characteristics: 25 cm? electrode surface; NAFION 
112 membrane; electrodes of coal cloth; and double serpentine 
topology in anode and cathode. Initial experimental conditions: 
26° C and 60% relative humidity. Experimental conditions: cell 
temperature is not controlled. 


2.2. Experimental procedure 


The PEMFC is fed with hydrogen and oxygen at the same 
pressure: 0.5 bar. The anode and cathode exhausts are kept 
closed. An electric load of 0.066 Q is connected. Periodically, 
CI are performed and an oscilloscope (Tektronix THS720) is 
used for cell voltage reading (2500 readings per CI test), with a 
sample period of 2 ms. 

In order to study the effect of the cathode flooding on the cell 
voltage and current, a data-acquisition card (PCLab PCL812G) 
and LabVIEW are used to record the cell voltage, with a sample 
period of 1s. This experimental data is used to estimate the 
amount of water produced in the electrochemical reaction (as it 
is discussed in Section 3.8). 


2.3. Uncertainty analysis 


Tektronix THS720 oscilloscope: dc voltage range, 400.0 mV 
to 880 V; de volts accuracy, + (0.5% of reading + 5 counts). 
e PCLab PCL&12G card: its measurement accuracy is 
+0.015% of the complete range. 

Electric load: it was calibrated using the HP 4192A 
impedance analyzer, which has a measurement error of 
+0.1%. 


3. Analysis and modelling 
3.1. Main assumptions 


The voltage drop across the fuel cell (Veen) can be written 
as a function of the steady-state open-circuit voltage of the cell 
(Eoc), and the over-voltages of the anode (na), the cathode (ne) 
and the membrane (7m) [21,26]: 


Vee = Eoc — Na — Nm — Ne (1) 


The open-circuit voltage (Eoc) is modelled as an ideal voltage 
source. The external operating conditions are kept constant dur- 
ing the experiments. Therefore, the over-voltages (Na, Nnm, Nc) 
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Fig. 1. PEMFC electric models: (a) complete model and (b) simplified model. 


are assumed to be only dependent of the FC internal phenom- 
ena. The over-voltages of the different layers can be modelled 
as electric components that represent the FC dynamic behavior. 


3.2. Equivalent circuit for a PEMFC 


The cell voltage in Eq. (1) can be modelled by the electric 
circuit shown in Fig. la, which is composed of Randles mod- 
els connected in series [22,18]. Zęķę and Zy are the Warburg 
impedances associated to the gas diffusion in the anode and the 
cathode, respectively. R5 and Rọ are the charge transfer resis- 
tances in the anode and the cathode. C4, and C4, are the double 
layer capacities in the anode and the cathode. Finally, Rm is the 
membrane resistance. 

Two additional hypotheses are made in order to simplify the 
cell model shown in Fig. la. As a result, the simplified model 
shown in Fig. 1b is obtained. 


(1) The oxygen reduction reaction in the cathode in very slow in 
comparison with the hydrogen oxidation reaction [27,28]. 
Therefore, the anode over-voltage is very small in compar- 
ison with the cathode over-voltage. As a consequence, the 
anode over-voltage contribution to the cell voltage can be 
neglected in the model. 

(2) The double layer capacity, which is usually represented by 
constant-phase elements, is represented by a pure, single- 
frequency theoretical capacity [19,16,22,29]. 


3.3. Modelling of the Warburg impedance 


The Warburg impedance (Zw) can be written in the Laplace 
domain as a function of the finite length diffusion [18]: 


tanh ./stq 
a/ STq 


where the diffusion resistance (Ra) and the diffusion time con- 
stant (tq) can be calculated from the following expressions: 


RTS 


Zw(s) = Ra (2) 


Ra = SST 3 

i SC Dn? F? a 
82 

= — 4 

Td D ( ) 


The following approximation of the Warburg impedance in 
Eq. (2) is proposed: 


Ri R3 


Zw(s) = 
WO) = TER Cs | TE RCs 


(5) 


Rn WE ye 
VW 


Fig. 2. PEMFC electric model with Warburg impedance calculated from Eq. 
(5). 


Table 1 

Fitted values of the parameters in Eq. (8) 

Parameter Value 

Ry (Q) 0.8463 

Rp (2) 0.1033 

Cı (F) 0.3550 

C2 (F) 0.03145 

where 

R* = R;Rą for i= {1,2} (6) 
x Td : 

Ci=C;— for i= {1,2} (7) 

Ra 


or, equivalently 


B R Ra 
Zw(s) = Ra + (8) 


14+ RiCitas 1+ RoCotas 


This approximation is equivalent to model the cell by using 
the circuit shown in Fig. 2. Four new parameters have been 
introduced in Eq. (8): R,, R2,C, and C2. The values of 
these four parameters are calculated by fitting Eq. (8) to Eq. 
(2) in the frequency range from 2 x 107! s7! to 1.5 x 10° 
s71, considering tq = 1. The calculated values are shown in 
Table 1. 

The exact value of the impedance, calculated from Eq. (2), 
and the approximated value, obtained from Eq. (8) using the 
previously calculated values of R1, R2, Cı and C2, are plotted 
in Fig. 3. Also, the absolute error is shown in Fig. 3. 


Joug 


10 10° 10 


Jor; 


Fig. 3. Warburg impedance calculated from Eq. (2) (—) and Eq. (8) (x). Abso- 
lute error. 
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Table 2 

Estimated parameters 

Parameter Description Unit 
Cai Double layer capacitance F 
Ra Diffusion resistance Q 
Rp Charge transfer resistance Q 

Ta Diffusion-related time constant s 
Rm Membrane resistance Q 


3.4. Modelling of the PEMFC impedance 


The cell impedance (Zeen), calculated from the circuit shown 
in Fig. 2, is the following: 


as? ++ bs +c 
Zeetl(S) = 7 (9) 


s + es? + fs+g 


where the parameters a, b, c, d, e, f and g can be calculated from 
the following expressions: 


a = 9.76 x 10-*Ryty (10) 
b = 0.304Rpta + 3.38 x 107° Rata (11) 
c = Rp + 0.949 Ra (12) 
d = 9.76 x 10 Ca Rpt (13) 


e = 3.37 x 107°Ca Rata + 0.3048Ca Rpta +9.76 x 10-475 


a4 
f= Cal Rp + 0.3048tg + 0.949Cq Ra (15) 
sad (16) 


The cell impedance described by Eq. (9) depends on the follow- 
ing four physical parameters: Car, Ra, Rp and ta. As it will be 
discussed later, the value of these parameters can be estimated 
from the cell response after CI. The estimated parameters are 
listed in Table 2. 


3.5. PEMFC response after CI 


The simplified model of the Warburg impedance proposed in 
Section 3.3 has been used (in Section 3.4) to obtain an analytical 
expression of the cell impedance. Now, this analytical expres- 
sion of the cell impedance will be used to obtain an analytical 
expression for the PEMFC voltage response after CI. 

The equivalent circuit of the CI experimental set-up is shown 
in Fig. 4a. The load impedance is Rjgaq. The cell impedance 


(a) Veen 


-4 ~ Fuel Cell 
N [eal % (b) E pe 
ameg AMMAM aw w 

‘ } R R Z 


Rioad ` —— 


Fig. 4. CI experimental set-up: (a) equivalent circuit and (b) dc impedance (Zr). 


(Zeen) is calculated from Eq. (9). As the cell impedance has 
only real poles, the cell response after CI is equal to the cell 
response after a current impulse at t = 0 (see Section 3.7): 


I(t) = slt) (17) 


where ô(t) is the Dirac delta, and Jo is the steady-state value of 

the current at tf = 07. It is calculated as follows: 
Veen (t = 07 

— Voenl ) (18) 

Zy(t = 07) 


where Zr is the dc impedance of the circuit shown in Fig. 4a. 
The equivalent circuit of the impedance Zr is represented in 
Fig. 4b. 

The cell voltage can be written, in the Laplace domain, as 
follows: 


Veels) = Voc — Zeen(s)L(s) (19) 


The time-domain expression of the cell impedance proposed 
in Section 3.4 is the following: 


3 2 
7 ar; + br; +c =f 
L [Zeen] = Rmd(t) + Y | H e 20 
[Zcen(s)] mô(t) = (ees eo) 


where r1, r2 and r3 are the three roots of the polynomial ds? + 
es + fs+g. 

The dynamic response of the cell voltage after CI is the 
following: 


Veen (t) = Voc — Io 
(Re 1 a /20? +a] eit +a et +03 a) 
ov 20 


(21) 
where 
ar? + br; +c 
7 3dr? + 2er, + f 


for i = 1,2,3 (22) 


Qi 
and also, it must be satisfied: 


CO 
I tg PP i= (23) 
-œ Oy 27 
The parameter o can be estimated from the CI experimental data 
as follows. The CI experimental data allow estimating the over- 
voltage due to the membrane resistance [30,21]. It corresponds 
to the abrupt rising of the voltage, from Vo to Vs, at time t = 0* 
(see Fig. 5). Vo is the initial voltage, and V, is the voltage at 
t = 0°. Both can be estimated from the CI experimental data. 
Then, the membrane resistance (Rm) can be calculated from the 
following expression: 


Ve — Vo = ln Rm (24) 


Finally, the dynamic response of the cell voltage can be cal- 
culated, for t > Ot (i.e., for Vee > Vs), from the following 
expression: 


Veeil(t) = Voc — Vs — Iola e" + age’ +a3e 3") (25) 
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Fig. 5. Estimation of over-voltage due to membrane resistance. 


where Voc is the steady-state, open-circuit voltage of the PEMFC 
(see Fig. 5). 

In order to illustrate the dynamic behavior of the simplified 
cell model shown in Fig. 2, the cell voltage after CI (calculated 
from Eq. (25)) and the impedance spectra are plotted in Fig. 6 
for selected values of the parameters Car, Ra, Rp and Ta. 


3.6. Model parameter estimation from CI data 


The voltage after CI of the cell represented by Eq. (9) is 
described by Eq. (25). The a; coefficients depend on the roots 
ri and on the parameters a, b, c, d, e and f (see Eq. (22)). These 
parameters depend on the physical parameters Car, Ra, Rp and 
Tq (see Eqs. (10)-(15)). As a consequence, Vee1(t) is a function 
of r1, r2, r3, Cai, Ra, Rp and ta. 

The analytical expression of Voey(t) as a function of 
r1, r2, r3, Cai, Ra, Rp and Tq is obtained by symbolic manipula- 
tion of Eqs. (25), (22), and (10)—(15). Finally, in order to estimate 
the value of these seven parameters (i.e., r1, r2, r3, Cai, Ra, Rp 
and Tq), this analytical expression of Voey(t) can be fit- 
ted to the experimental data of the cell voltage after CI 
(for Veen > Vs). 
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Fig. 7. Input signals: (a) interruption step and (b) unit impulse. 


The procedure described in this section allows estimating the 
value of the parameters Rm, Cai, Ra, Rp and tg. These parame- 
ters completely define the simplified model of the fuel cell shown 
in Fig. 2. 


3.7. Interruption step and impulse responses 


Consider a plant transfer function described by Eq. (26). It 
will be demonstrated that the plant response to the following two 


input signals is the same: an interruption step and a unit impulse 
(see Fig. 7): 


Fs) = 3 where ki, q eRandk; ti >0 (26) 
F 1+ kts eas as 


The plant response to the interruption step can be calculated 
by subtracting the response to a unit step, u(t), from that to 
f(t) = 1 Ge., the interruption step can be written as 1 — u(f)). 

The plant response to a unit step is the following: 


=j 1 _ =1 _ —t/ki Ti 
[Foi] =g DE + kts) a2 cage 
(27) 


The plant response to f(t) = 1 is equal to the plant steady- 
state response to a unit step: 


Jim Et [ro] a = lim 1 {SFOs \ = lim{F(s)} = ki 


(28) 
©) h 
0.996 | 
= 0992} if 
D | 
oD 
& 0.988) 
[s] |! 
> 
0.984 | 
E Re(Z X9) 
0.98 -# 
0 0.05 0.1 0.15 0.2 
time (s) 


Re(ZX9) 
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Fig. 6. Simulated cell voltage after CI and impedance spectra for selected parameter values. Baseline: tq = 0.5 s, Ca = 0.1 F, Ra = 6 mQ and Rp = 4 mQ. (a) 
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Fig. 8. Overall experiment results: (a) cell voltage and (b) cell current and estimated amount of water. 


Therefore, the plant response to an interruption step can be initial amount water in the cell: 


calculated by substracting Eq. (27) from Eq. (28): m 


Si — Poud =e) = ky e/t (29) (ALON = (xmo) + È (Axmo); (33) 


i=1 


On the other hand, the plant response to a unit impulse is the 4. Results and discussion 

following: 

Cathode flooding and drying have a relevant effect on the cell 

z5 k; ent it (30) performance. In order to study the effect of the cathode flooding 

1+ k Tis 7 on the cell model parameters, the parameter estimation method- 

ology proposed in Section 3.6 can be repeatedly applied during 

As Eq. (30) is equal to Eq. (29), the plant described by Eq. the flooding process. The experimental set-up was described 

(26) has the same response to an interruption step than to a unit in Section 2. The data obtained from each CI experiment run, 


THF) = L~! D+ kas 


impulse. corresponding to a given cathode flooding level, is used to 
estimate the parameters of the cell electric model shown in 
Fig. 2. 


3.8. Amount of water in the fuel cell The time evolution of the cell voltage and current are plotted 


in Fig. 8a and b, respectively. The voltage and current decrease 
as time passes, due to the cathode flooding. 

As the cathode and the anode exhausts are closed, all the water 
produced by the electrochemical reaction is accumulated inside 
the cell. This can be used to estimate the water accumulated 
inside the cell (xy,0). The initial mass of liquid water inside 


The amount of water (AxH,0) generated by electrochemi- 
cal reaction in the fuel cell during the time period Tm can be 
calculated from Eq. (31) [31]. The cell current (Zins), which 
is supposed constant during that time period, can be calcu- 
lated from the cell voltage and the load resistance (see Eq. 


(32)): the cell is zero. Also, the initial cell temperature is low enough 
C Meats, to consider that the amount of water in gas phase is negligible. 
Axm0 = S$ Nan "S (31) The method used to estimate the water generated in the cell 
ae was discussed in Section 3.8. The time evolution of the water 

Vins accumulated inside the cell is plotted in Fig. 8b. 
fins = Rioad (32) Several CI runs are performed for different flooding condi- 
tions. The measured and fitted cell voltages for the 8th, 10th 
After N time intervals, the total amount of water in the cell, and 12th CI runs are shown in Fig. 9a. The corresponding 


(XH20) y> can be calculated from Eq. (33), where (XH,0)p is the circuit models are calculated by applying the proposed method- 


~ 
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Fig. 9. Results from selected CI runs: (a) measured and fitted cell voltages obtained from 8th, 10th and 12th CI runs; and (b) impedance spectra calculated from 8th 
(o), 10th (x) and 12th (e) CI runs. 
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Table 3 
PEMFC model parameters calculated for selected CI runs 


CI run Ca (F) Ra (m2) Rp (mQ) Tq (S) Rm (mQ) Xm0 (g cm~?) SSE 
3rd 1.006 4.63 0.50 0.271 14.6 0.037 0.68 x107? 
7th 0.998 5.21 1.59 0.406 13.8 0.112 0.72 x107? 
8th 0.994 5.62 2.01 0.463 13.9 0.120 1.01 x107? 
9th 0.991 5.80 3.88 0.551 10.4 0.129 1.78 x107? 
10th 0.989 6.35 5.30 0.779 13.9 0.141 1.03 x107? 
11th 0.986 8.47 5.54 0.821 12.4 0.150 1.63 x107? 
12th 0.984 10.8 5.71 0.874 14.6 0.167 1.06 x107? 
13th 0.973 18.0 6.20 0.947 17.7 0.187 1.51 x107? 
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ology, and they are used to obtain the impedance spectra, which 
are plotted in Fig. 9b. As expected [16,32], the width of the 
impedance arcs increases as a result of the cell flooding. 

The model parameters calculated from the data of some CI 
runs, using the Matlab curve fitting toolbox, are shown in Table 3. 
The parameter SSE is the sum of the quadratic error of the curve 
fitting. 

Note that tg, Ra, Rp and Cay depend on the amount of water 
accumulated inside the cell (7,0). In particular, the diffusion 
resistance (Ra) depends on the diffusion media coefficient and 
the oxygen concentration on the active surface (see Eq. (3)). 
The relationship between Ra and xy,0 is plotted in Fig.. The 
points calculated from the 7th, 8th, 9th, 10th, 11th, 12th and 
13th CI runs are represented by the symbol x. The third-order 
polynomial fitted to these points is the following: 


Ra = 25Xi1,0 — 8-3Xfpo + 0.96xmo — 0.033 (34) 


5. Conclusions 


A PEMFC circuit model, well suited for control applications, 
and a methodology for model parameter estimation from CI 
data, have been proposed. A key advantage of this approach 
is that the equipment required to perform the CI tests is easily 
portable and inexpensive. As a consequence, it can be imple- 
mented in cell commercial systems for in-field assessment of 
the cell. The proposed modelling methodology has been suc- 
cessfully applied to the study of the cathode flooding pheno- 
menon. 
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